A New Feature in Some Quasi-discontinuous Systems * 



WU Shun-Guang( ) 1 - 2 HE Da-Ren( ) 1 

1 Complexity Science Center, Yangzhou University, Yangzhou 225002 
2 Institute of Low Energy Nuclear Physics, Beijing Normal University, Beijing 100875 

(Received 20 October, 1999) 



ON 
ON 
G\ 

> 

o 

(N 



o : 
o ■ 

(N ■ 

on : 

ON . 

X3 ■ 
i 1 

o : 

^ ■ 

O ■ 



Many systems can display a very short, rapid changing 
stage (quasi-discontinuous region) inside a relatively very long 
and slowly changing process. A quantitative definition for the 
" quasi-discontinuity" in these systems has been introduced. 
We have shown by a simplified model that extra-large Feigen- 
baum constants can be found inside some period-doubling 
cascades due to the quasi-discontinuity. As an example, this 
phenomenon has also been observed in Rose-Hindmash model 
describing neuron activities. 
PACS: 05.45. +b 

Recently, there has been considerable interest in piece- 
wise smooth systems (PWSSs). Such models usually de- 
scribe systems displaying sudden, discontinuous changes, 
or jumping transitions after a long, gradually varying 
process. These systems may show some behaviors appar- 
ently different from those of the everywhere - differen- 
tiable systems (EDSs) [1-5]. In fact, the sudden changes 
in the above processes also need time. Therefore, such a 
process can be everywhere smooth if one describes it with 
a high enough resolution. Usually, in the largest part of 
the process, a quantity changes very slowly. It has a dras- 
tic changing only in one or several very small stages. We 
suggest to call the stage as a " quasi-discontinuous region 
(QDR)" and shall define a "quasi-discontinuity (QD)" 
inside it quantitatively. A system that can display QDR 
in its processes may be called a " quasi-discontinuous sys- 
tem (QDS)" . Obviously, QDS is a much wider conception 
than PWSS and may serve as an intermediate between 
EDS and PWSS. 

In order to show our basic idea and the first character- 
istic of QDS, we have constructed a model map as shown 
in Fig.l. The map reads: 



f hi* 

fix) = < f 3 (x 
{ h(x 



ki(x - xi) + yi 
A (x-x ) 2 + y 
^/r 2 - (x - o x ) 2 + o y 
k 2 (x - Xi) + j/4 
k 3 (x - 1) 



x e [o,x{), 
x e [2:1,2:3), 

X G [2:3,2:4), 
x G [2:4,2:5), 

x e [x 5 , 1]. 



(1) 



As can be seen in Fig.l, the slope of /1 branch is a unit. 
It is simulating the slowly changing part of the process. 
Branch fi is a linear line with a very large negative slope 
k 2 . Branch is a small part of a circle introduced for 
a smooth connection of f 2 and fi. The center of the 
circle locates at (o x ,o y ), and its radius is r. Branches 



/3 and /4 can simulate the small drastic changing part. 
In Eqs. (1), A is chosen as the control parameter. It 
is obvious that the fixed point at P 2 undergoes period- 
doubling bifurcation when A changes inside a certain pa- 
rameter range. (xj,yj) denote the coordinates of points 
Pj (j=0,...,5), respectively. They are determined by the 
conditions of smooth connections between neighboring 
branches. For certain function forms of f 2 and fi, the 
circle of still may be very large or small. We define 
another parameter a to fix it. Therefore, o x , o y , r, k, 
and (xj,yj) are all functions of A and a. Their explicit 
forms will not be shown in this short letter. The param- 
eter ranges chosen for this study are A G [8.0,9.0] and 
a G [0.99995,0.999997]. 

Now we will define QDR and QD in this model. Ac- 
cording to the geometrical properties of Eqs. (1), one can 
obtain the following conclusions. When a = 1, r = 0, the 
first order derivative of the map function is discontinu- 
ous at 2:3 = 2:4. The second order derivative shows a 
singularity, that is, an infinitely large value here. When 
a 6 (0, 1), the branch fa has a finite length. The first or- 
der derivative of the map function is continuous at both 
X3 and X4. The second order derivative value between 
them is finite but trends toward infinite when a — > 1. In 
this case, the maximum value of the second order deriva- 
tive of the map function between X3 and 2:4 may be used 
to describe the "quasi-discontinuity (QD)". So we shall 
define QD as 
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and define QDR as 



A =11(2) -a:(l)|, (3) 
where a:(2) and 2;(1) are between £3 and 2:4, and satisfy 
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Acording to the definations (2) and (3), the QDR and 
QD for Eqs. (1) can be expressed as 



x(l) = o x 



kor 



V2 + fc 



x(2) = Xi, 



and 



r2 — (xi — o x ) 2 ] 3 / 2 ' 
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respectively. 

When a = 1 it is reasonable to observe one of the 
typical behaviors of PWSS. That is the interruption of a 
period - doubling bifurcation cascades by a type V inter- 
mittency [2,3]. 

When a is smaller than, but close to 1, there is a 
QDR between x% and X4 instead of the non-diffcrentiablc 
point. The mapping is everywhere smooth, so the period- 
doubling bifurcation cascade should continue to the end. 
However, there is a drastic transition of the mapping 
function slope in a very small QDR that makes all fur- 
ther bifurcation points compressed into a relatively much 
shorter parameter distance. The Feigenbaum constants 
Si (i=3,4,5 or even more), influenced by the compression, 
should show some extraordinary values. That is exactly 
what we have observed. Table 1 shows the data about 
three cascades. In the table, n indicate the sequence 
number of doubling, S n (i) (i = 1 — 3) are the Feigenbaum 
constants of the cascade i. The parameter values a(i) and 
the maximum value of the QD, for each cascade are 
indicated in the caption. In the table 5 no data are ob- 
tained from Ref. [6]. They are listed here for a compar- 
ison with the corresponding ones obtained in a typical 
everywhere smooth situation. As can be seen in table I, 
when is large, a lot of Feigenbaum constants, £3, £4, 
£5, <5g and Si are extraordinary. The further constants 
may be considered as ordinary, but they converge to the 
universal Feigenbaum number very slowly. When is 
smaller, only £3, 64 and S$ are apparently extraordinary. 
The further constants converge much faster. When 
is very small, the whole Feigenbaum constant sequence is 
very close to the standard S no data. That may indicate a 
smooth transition from QDS to a EDS. Also, from these 
data one can believe that the extraordinary Feigenbaum 
constants in the period-doubling cascades are induced by 
QD of the system. Based on this understanding we sug- 
gest the use of the common extraordinary Feigenbaum 
constant 5j(j = 3,4) to signify this phenomenon. The 
relationship between k, the QD, and the symbol of the 
phenomenon Sj (j = 3,4), have been computed. Figure 
2 shows the result of function k — Sj ( Although k is de- 
pendent on both a and A, our numerical results demon- 
strate that k is not sensitive to the parameter A at a 
given a. Therefore, it is possible to choose the maxmum 
k to represent QD of a whole diagram. For example, for 
the bifurcation points A n (n — 0, 1, 2, 12), indicated in 
the second column of Table 1, the corresponding k are 
220.9 x 10 6 , 208.2 x 10 6 , 205.6 x 10 6 , 205.0 x 10 6 , 
respectively. So we choose 220.9 x 10 6 as the represen- 
tative k of the bifurcation diagram). One can see that 
S4 increases, but 63 decreases when k becomes larger and 
larger. 

It is important to find examples of this kind of interest- 
ing phenomenon in practical systems. We have done such 
a study in Rose-Hindmarsh (R-H) model. The model, 
which describes neuronal bursting [7], can be expressed 



by 

= y - ax 3 + bx 2 + I - z, 
% = c-dx 2 -y, (6) 
f t =r[s(x-x*)-z], 

where x is the electrical potential of the biology mem- 
brane, y is the recovering variable, z is the adjusting 
current, a, b, c, d, s and x* are constants, r and / are cho- 
sen as the control parameters. We shall take a — 1, b = 
3, c = l,d= 5,x* = 1.6, s = 4 for this study. Fig. 3 shows 
the Poincare map of a strange attractor observed when 
/ = 2.9, r = 0.00433 (Here the Poincare section is de- 
fined as the coordinate value of z axis at the maximum 
in x direction of the trajectory. We have also tested some 
different definitions of Poincare section, the results have 
shown that all of them are qualitatively the same as each 
other) . It is clear that the iterations in the region \z\ , Z2] 
change very rapidly. Therefore, we call this region as a 
QDR. 

Table 2 shows the critical bifurcation parameter val- 
ues and the corresponding Feigenbaum constants for a 
period-doubling bifurcation cascade. One can see that Si 
and 62 are larger than ordinary values. As our compu- 
tation has confirmed, that means an interruption of the 
cascade by a collision of the periodic orbit with the QDR 
when the first time period-doubling finished. 

For a comparison with the function k — Si shown in Fig. 
2, we have computed the bifurcation diagrams with / = 
2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, and r e [0.1 x 10~ 2 , 4x 10" 2 ]. 
The results are shown in Fig. 4 (Here, we also choose the 
maxmum k to represent the quasi-discontinuouty of one 
bifurcation diagram). They are in a qualitative agree- 
ment with those in Fig. 2. 

In conclusion, we have found some extraordinary 
Feigenbaum constants in some period-doubling bifurca- 
tion cascades in a constructive and a practical system. 
The mechanism of the phenomenon is that a periodic or- 
bit near a critical point of bifurcation crosses a QDR in 
the system. This understanding may be important for 
the experimental scientists because very often they can 
measure only the first several Feigenbaum constants in 
a real experiment. After observing strange Feigenbaum 
constants, they can verify if their system is a QDS with 
the knowledge in this discussion. Moreover, our results 
also demonstrate that between typical PWSSs and EDSs 
there can be a type of transitive systems. 
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TABLE I. The data for period-doubling bifurcation cas- 
cades in Eqs. (1) with a = 0.9999970 and n max = 220.9 x 10 6 
for cascade (1); a = 0.9999196 and K max = 8.2 x 10 6 for cas- 
cade (2); a = 0.9987537 and n max = 0.5 x 10 6 for cascade (3). 
The data are obtained by numerical method with a quadruple 
precision. 
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Mi) 


M2) 


M3) 


Sn 


1 


4.626168416 


4.626168416 


4.626168416 


4.744309468 


2 


4.638635250 


4.638635250 


5.965311262 


4.674447827 


3 


8.823236594 


9.982064249 


7.746315906 


4.670792250 


4 


16.553602405 


11.112580986 


4.756159889 


4.669461648 


5 


1.619920263 


3.966760226 


4.705758197 


4.669265809 


6 


4.147211067 


4.411839476 


4.674983020 


4.669214270 


7 


5.173386810 


4.526007508 


4.670687838 


4.669204451 


8 


4.765263421 


4.709611842 


4.669607200 


4.669202201 


9 


4.757647408 


4.671017586 


4.669986034 


4.669201737 


10 


4.687625667 


4.677922775 


4.669375410 


4.669201636 


11 


4.701929236 


4.703838840 


4.699239861 


4.669201614 
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FIG. 1. Function of Eqs. (1) when A = 8.75407, a = 0.998. 
The inset is the enlargement area near P3 and P4,. 




TABLE II. The data of the bifurcation critical parameter 
values and the corresponding Feigenbaum constants observed 
in system (6) with I = 2.8. 
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5n 





3.52000 
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1.55550 


19.84344 
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1.45650 


6.14905 


3 


1.44040 


5.19350 


4 


1.43730 


4.76981 


5 


1.43665 
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FIG. 2. The function k — Si for Eqs. (1) with the param- 
eter value a G [0.99966, 0.999997]. 
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FIG. 3. The Poincare map of R-H model. Here the range 
[zi, Z2] schematically indicates the QDR, the real QDR calcu- 
lated via Eq.(3) is too small to be clearly expressed. 
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FIG. 4. The function k — Si for period-doubling bifurca- 
tion's in R-H model. 
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